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Abstract 

In this paper we prove two results regarding reconstruction from magnitudes of frame coefficients (the 
so called ’’phase retrieval problem”). First we show that phase retrievability as an algebraic property 
implies that nonlinear maps are bi-Lipschitz with respect to appropriate metrics on the quotient space. 
Second we prove that reconstruction can be performed using Lipschitz continuous maps. Specifically we 
show that when nonlinear analysis maps a, j3 ■. H ^ R"* are injective, with a{x) — {\{x, fk)\)T=i s-nd 
l3{x) = {\{x, fk)\^)T=ij where {/i,...,/m} is a frame for a Hilbert space H and H = H/T^, then a is 
bi-Lipschitz with respect to the class of ’’natural metrics” Dp{x, y) = miriip ||® — whereas (1 is bi- 

Lipschitz with respect to the class of matrix-norm induced metrics dp(x, y) = IRx* — j/i/* ||p. Furthermore, 
there exist left inverse maps ujjip ■ R"* H of a and /3 respectively, that are Lipschitz continuous with 
respect to the appropriate metric. Additionally we obtain the Lipschitz constants of these inverse maps 
in terms of the lower Lipschitz constants of a and /3. Surprisingly the increase in Lipschitz constant is a 
relatively small factor, independent of the space dimension or the frame redundancy. 


1 Introduction 

Assume T = {/i, / 2 , • ■ ■ j /m} is a frame (that is a spanning set) for the n-dimensional Hilbert space H. In 
this paper H can be a real or complex Hilbert space. The results in Section [3] apply to both cases, and the 
constants have the same form. 

On H we consider the equivalency relation a; ~ y iff there is a scalar a of magnitude one, |a| = 1, so that 
y = ax. Let H = H/ ^ denote the set of equivalence classes. Note H \ {0} is equivalent to the cross-product 
between a real or complex projective space of dimension n — 1 and the positive semiaxis R"*". 

In this paper we use x to denote the equivalency class of x in H . Nevertheless, for simplicity, x is used 
in place of x when there is no ambiguity. 

Let a and /3 denote the nonlinear maps 


a-. H 

, «(*) = (|(a;,/fc)|)i<fc<^. 

(1) 

/3 : 7L R”" , 

/3(x) = 

(2) 


1 


The phase retrieval problem, or the phaseless reconstruction problem, refers to analyzing when a (or /3) 
is an injective map, and in this case to finding ’’good” left inverses. 

The frame is said to be phase retrievable if the nonlinear map a (or /3) is injective. In this paper 
we assume a and /3 are injective maps (hence T is phase retrievable). The problem is to analyze Lipschitz 
properties of these nonlinear maps, and then to extend the unique left inverse from the image of H through 
the nonlinear maps a, /?, to the entire space K.™ so that they remains Lipschitz continuous. 

A continuous map / : {X, dx) (Y, dy), defined between metric spaces X and Y with distances dx and 
dy respectively, is Lipschitz continuous with Lipschitz constant Lip(/) if 


Lip(/) := 


sup 

Xi ,X2€X 


dy{f{xi),f{x 2 )) 

dxixi,X2) 


< oo. 


The map / is called bi-Lipschitz with lower Lipschitz constant a and upper Lipschitz constant b if for every 
Xi,X2 £ X, 

adx{xi,X 2 ) < dy{f{xi),f{x 2 )) < bdx{xi,X 2 ). 

Obviously the smallest upper Lischitz constant is 6 = Lip(f). If / is bi-Lipschitz then / is injective. 

The space H admits two classes of inequivalent metrics. We introduce and study them in detail in section 
[21 In particular consider the following two distances: 


D2{x,y) 

= mm 2 = + Ibf - 2|(x, y)| 

(3) 

di{x,y) 

= \\xx*-yy*\\, = J{\\xf + \\y\\y-4\{x,y)\^ 

(4) 


When the frame is phase retrievable the nonlinear maps a : (iL, D 2 ) —>■ IMI2) ■ {dd,di) — >■ 

(R™, 11 - 112 ) are shown to be bi-Lipschitz. This statement was previously know for the map (3 in the real and 
complex case (see 0121 ITTjl. and for the map a in the real case only (see [24l [11] [13]). In this paper we 
prove this statement for a in the complex case. Denote by Oq and ap the lower Lipschitz constants of a and 
P respectively. In this paper we prove also that there exist two Lipschitz continuous maps w : (M™, 11-112) —^ 
{H,D 2 ) and ip : (R™, II-II 2 ) — {H,d{) so that uj{a{x)) = x and ip{P{x)) = x for every x G H. Furthermore 
the upper Lipschitz constants of these maps obey Lip(a;) < ^^5 Lip('0) < Surprisingly this shows 

the Lipschitz constant of these left inverses are just a small factor larger than the minimal Lipschitz constants. 
Furthermore this factor is independent of dimension n or number of frame vectors m. 

The organization of the paper is as follows. Section [2] introduces notations and presents the results for 
bi-Lipschitz properties. Section [3] presents the results for the extension of the left inverse. Section [4] contains 
the proof of these results. 


2 Notations and Bi-Lipschitz Properties 

On the space H we consider two classes of metrics (distances) induced by corresponding distances on H and 
respectively: 

1. The class of natural metrics. For every 1 < p < 00 and x,y £ H define 

Dp{x,y) = min \\x - ay\\ (5) 

|a| = l 

When no subscript is used, ||-|| denotes the Euclidian norm, ||-|| = ||-|| 2 . 
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2. The class of matrix norm induced metrics. For every 1 < p < oo and x,y € H define 


dp{x,y) 


\\xx* 



for 1<P<00 
maxi<fc<„ (Tfc for p = oo 


(6) 


where {crk)i<k<n are the singular values of the operator xx* — yy*, which is of rank at most 2. 

Our choice in ([51) corresponds to the class of Schatten norms that extend to ideals of compact operators. 
In particular p = oo corresponds to the operator norm IHlgp in Sym(iJ) = {T : H ^ H , T = T*}; p = 2 
corresponds to the Frobenius norm ll’llii’r in Sym(iJ); p = 1 corresponds to the nuclear norm H-H^ in Sym(iJ): 

doo{x,y) = \\xx*-yy*\\^p , d 2 {x,y) = \\xx* - yy*\\p^ , di{x,y) = \\xx* - yy*\\.. 

Note the Frobenius norm ||T||p^ = •\/trace(TT*) induces an Euclidian metric on Sym(iJ). In [7] Lemma 
3.7 we computed explicitly the eigenvalues of |x, t/]. Based on these values, we can easily derive explicit 
expressions for these distances: 

d^{x,y) = l|xf - \\yf I + l\/{\\xf + \\y\\y-i\{x,y)\^ 

d 2 {x,y) = ^J\\x\f + \\y\f -2\{x,y)\^ 

diix,y) = \J{\\xf + \\y\\^Y - A\{x,y)\^ 

Since xx* — yy* = Ju,u] for a; = + v) and y = ^{u — v). 

To analyze the bi-Lipschitz properties, we define the following three types of Lipschitz bounds for a. 
Note that the Lipschitz constants are square-roots of those constants. 


(i) The global lower and upper Lipschitz bounds, respectively: 


^0 = inf 

x,y£H 


\\oi{x) — a 


D2{x,yY 


Bq = sup 

x,y£H 


||Q:(a;) - a{y)\\l ^ 
D2{x,yY 


(ii) The type I local lower and upper Lipschitz bounds at z € H, respectively: 


A(z) = lim inf 

’ ^0 x,yeH 

D2{x,z)<r 
D2{y,z)<r 


||a(a:) - aiy)\\l 
D2{x,yY 


B{z) = sup inf _ 

r—vO x,yCiH 
D2{x,z)<r 
D2{y,z:)<r 


||a(a:) - aiy)\\l _ 
D2{x,yY 


(hi) The type II local lower and upper Lipschitz bounds at z G H, respectively: 


A(z) = lim inf 

'■->■0 x€H 

D2(x,z)<r 


||a(x) - a{z)\\l 

D2{x,zY 


B{z) = sup inf 

r^O x^H 

D2{x,z)<r 


\\a{x) - a{z)\\l 
D2{x,yY 


Similarly, we define the three types of Lipschitz constants for j3. 
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(i) The global lower and upper Lipschitz bounds, respectively: 


Oq 


inf 

x,y£H 


di{x,yy 


bo = sup 

x,y£H 


mx)-/ 3iy)r, 

di{x,y)‘^ 


(ii) The type I local lower and upper Lipschitz bounds at z € H, respectively: 


a{z) = lim inf 

'’->■0 x , y&H 
di{x,z)<ir 
dl{y,z)<r 


W{x)-P{y)\\l 

di{x,yY 


b{z) = liin 

r^O 


sup 

x,y£H 
di {x,z)<ir 
diiv,z)<r 


mx)-l3iy)\\l 

di{x,yY 


(hi) The type II local lower and upper Lipschitz bounds at z € H, respectively: 


a(z) = lim inf 

’’->■0 y€H 

di{x,z)<.r 


mx)-(3iz)\\l 

dl (x, z)2 


b(z) = lim 

^ ^ r-)-0 


sup 

y^H 

di{x,z)<.r 


WPjx) - I3{z)\\l 
di{x,zY 


Note that due to homogeneity we have Aq = A(0), Bq = B{0), qq = a(0), bo = 6(0). Also, for z ^ 0, we 
have A(z) = A(z/ ||zl|), B{z) = B{z/ || 2 ;||), a{z) = a^zj || 2 ||), b{z) = b{z/ \\z\\). 

We analyze the bi-Lipschitz properties of a and /3 by studying these constants. 


2.1 Bi-Lipschitz Properties of a 

The real case H = R" is studied in m- We summarize the results as a theorem. 

Recall that B = {fi, ■ ■ ■ , /m} is a frame in H if there exist positive constants A and B for which 

m 

A\\xf<Y,\{x,h)\^<B\\xf. (7) 

k=l 

We say A [resp., B] is the optimal lower [resp., upper] frame bound if A [resp., B] is the largest [resp., 
smallest] positive number for which the inequality ([7|) is satisfied. 

For any index set / C {1, 2, • • ■ , m}, let iF[I] = {fk, k G 1} denote the frame subset indexed by I. Also, 
let <Ji[I] and a‘^[I] denote the upper and lower frame bound of set B[I], respectively. That is, 

CTi [/] = Amax ( ^ fkfk ) and al [/] = A^in ( XI 
\kel ) \kGl 

Theorem 2.1 (la) Let IF C K” be a phase retrievable frame for R”. Let A and B denote its optimal 
lower and upper frame bound, respetively. Then 

(i) For every t) ^ x G R", A{x) = a-'^{supp{a{x)) where supp{a{x)) = {k, (x,fk) fz 0}; 

(ii) For every x G R”, A = A; 

(hi) Ao = A(0) = min7c{i,2,...,m}(o-n[l] > 0; 
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(iv) For every x G K", B{x) = B{x) = B; 

(v) Bo = B{0) = B{0) = B. 


Now we consider the complex case H = C". We analyze the complex case by doing a realification first. 
Consider the R-linear map j : C" —>■ defined by 




real( 2 :) 

imag(z) 


This realification is studied in detail in [7]. We call j{z) the realification of z. For simplicity, in this paper 
we will denote ^ = j(x), r? = j(2/), C = j(2), = j(/), = j(d), respectively. 

For a frame set T = {/i, / 2 , • • • , fm}, define the symmetric operator 


where 

is a matrix in ]R 2 nx 2 n^ 

Also, define 5 : ^ Sym(R2") by 


S{0 


E 


1 




( 8 ) 


We have the following result: 

Theorem 2.2 Let T C C" he a phase retrievable frame for C". Let A and B denote its optimal lower and 
upper frame bound, respetively. For any z € C", let ( = j(z) be its realification. Then 

(i) For every 0 ^ z G C”, A(z) = A2n-i(5(C)) ; 

(ii) Ao = A(0) > 0 ; 

(hi) For every z G C", i(z) = \ 2 n -1 (s{C) + T,k-.{zM=o f 

(iv) i(0) = A ; 

(v) For every z G C", B{z) = B{z) = Ai ^5(0 + Y.k-.{zM=o > 

(vi) Bo = B{0) = B{0) = B . 
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2.2 Bi-Lipschitz Properties of f3 

The nonlinear map /3 naturally induces a linear map between the space Sym(i7) of symmetric operators on 
H and K™: 

.4 : Sym(it) ^ R™ , AiT) = ({ThJk))i<k<m (9) 

This linear map has first been observed in m and it has been exploited successfully in various paprs e.g. 

mulls]. 

Let SP’‘^{H) denote the set of symmetric operators that have at most p strictly positive eigenvalues and 
q strictly negative eigenvalues. In particular denotes the set of non-negative symmetric operators 

of rank at most one: 

= {xx*, x€H} (10) 

where x* : H ^ <C, x*{y) = {y,x) is the dual map associated to vector x. In [7] we studied in more depth 
geometric and analytic properties of this set. In particular note I5{x) = .4(|a;,a;]) where 

= + yA") (ii) 

denotes the symmetric outer product between vectors x and y. The map P is injective if and only if A 
restricted to is injective. 

In previous papers mill] we showed the following necessary and sufficient conditions for a frame to give 
phase retrieval. 

Theorem 2.3 (Bin]) The following are equivalent: 

(i) The frame T is phase retrievable; 

(ii) ker{A)^S^’^{H) = {{)}; 

(hi) There is a constant oq > 0 so that for every u,v £ H 

^ m 

(iv) There is a constant > 0 so that for every x^y ^ H, 

\\p{x)-P{y)f >ao{di{x,y) f (13) 

In |7] we also showed a theorem that can be restated as follows: 

Theorem 2.4 If T is phase retrievable, then there is a constant bo > 0 such that for every x,y € H, 

\\P{x) - P{y)f <bo{di{x,y))^ (14) 


||uf ||uf - (imag((u, ?;)))" 


( 12 ) 


where bo is given by 

m m 

= II Ii^if^i X (i'eal((a:,/fc)(/fc,2/)))^ = max ^ \{x,fk)f = ll'r||s(;2 ,4) 
k = y =1^ 

k—1 fc=l 

In the last expression T : H ^ C'” is the analysis operator defined by x ^ /fc))fcLi ■ 


(15) 
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Remark 2.5 An upper bound of bo is given by 


bo<B( max \\fk\\] < (16) 

yl<fe<?7i / 

where B is the upper frame bound of T. 

We give an expression of the local Lipshitz bounds as well. Define 7^ : —?► Sym(R^") by 

m 

k=\ 

Theorem 2.6 Let B be a phase retrievable frame for H = C”. For every 0 z £ H, let f = '^{z) denote 
the realification of z. Then 

(i) aiz) = d{z) = \2n-im))/Kf; 

(ii) b{z) =b{z) = \i{TZ{C))/ WCf; 

(hi) ([7]) a(0) = ao = min||<.||=i X2„-i (7^(C)); 

(iv) 5(0) = min||^||=i X;r=i \{xJk)fi 

(v) 6(0) = 6(0) = bo; 


3 Extension of the Inverse Map 

All metrics Dp and dp induce the same topology as shown in the following result. 

Proposition 3.1 (i) For each 1 < p < oo, Dp and dp are metrics (distances) on FI. 

(ii) (Dp)i<p<oo ore equivalent metrics, that is each Dp induces the same topology on FI as Di. Additionally, 
for every 1 < p,q < oo the embedding i : (H,Dp) —>■ [H,Dq), i{x) = x, is Lipschitz with Lipschitz 
constant 

= max(l,n5"p). (18) 

(iii) For 1 < p,q < oo, (dp)i<p<oo ore equivalent metrics, that is each dp induces the same topology on iJ 
as di- Additionally, for every 1 < p,q < oo the embedding i : {H,dp) —>■ {H,dq), i{x) = x, is Lipschitz 
with Lipschitz constant 

= max(l,2?”p). (19) 

(iv) The identity map i : {F[,Dp) —>■ [FI,dp), i[x) = x, is continuous with continuous inverse. However it 
is not Lipschitz, nor is its inverse. 
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(v) The metric space {H,Dp) is Lipschitz isomorphic to endowed with Schatten norm IHIp. The 

isomorphism is given by the map 

: A ^ S«(il) , K„(.) = { ^ ^ (20) 

0 ij X = 0 

i—1 

The embedding is bi-Lipschitz with the lower Lipschitz constant min (22 p,np =) and the upper 
Lipschitz constant •\/2ma.x{n^~p,2p~^). In particular for p = 2 the lower Lipschitz constant is 1 and 
the upper Lipschitz constant is \f2. 

(vi) The metric space (H,dp) is isometrically isomorphic to endowed with Schatten norm ll’Hp. 

The isomorphism is given by the map 

up : H ^ , Kp{x) = XX*. (21) 

In particular the metric space {H,di) is isometrically isomorphic to endowed with the nuclear 

norm 

(vii) The nonlinear map l : {H,Dp) {H,dp) defined by 

f —0= if X pi 0 

,(x) = \ VM (22) 

0 ij X = 0 

i_i i_i. 

is bi-Lipschitz with the lower Lipschitz constant min (22 p ^np and the upper Lipschitz constant 
v2max(n2 p,2p ^). 

Remark 3.2 (i) Note the Lipschitz bound Lp^^^ is egual to the operator norm of the identity between 

(C-, IMIp) and (C", INI,).- = ||/|bp(cp)Uzo(Cp)- 

(ii) Note the eguality Lp^^ ,^ = 

The results in Section [2l together with the previous proposition, show that frame T is phase retrievable 
then the nonlinear map ([T]) [resp., (H))] is bi-Lipschitz between metric spaces {H,Dp) [resp., {H,dp)] and 
(R™,||-||^). Recall that the Lipschitz constants between (iL, £> 2 ) [resp. (iL, di)] and (R’",||-|| = 11-112) are 
given by [resp., and [resp., Vbo]: 

^/A^D 2 {x, y) < ||Q((a;) - a{y)\\ < y^D 2 {x, y) (23) 

^/^dl{x,y) < \\l3{x)- P{y)\\ < ^/b^di{x,y) (24) 

Clearly the inverse map defined on the range of a [resp., /3] from metric space {a{H), H-H) [resp., {/3{H), H-H)] 
to {H,D 2 ) [resp., iH,di)]: 

Cj : a{H) C K"* —>■ II , <I'(c) = a; if a(x) = c (25) 

fj : fI{H) , fjic) = a; if fl{x) = c (26) 

is Lipschitz with Lipschitz constant [resp., In this paper we prove that both w and ip can be 

extended to the entire R™ as a Lipschitz map with Lipschitz constant that increases by a small factor. 

The precise statement is given in the following Theorem which is the main result of this paper. 
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Theorem 3.3 Let JF = {/i,..., fm} be a phase retrievable frame for the n dimensional Hilbert space H, 
and let a, P : H ^ R'” denote the injective nonlinear analysis map as defined in m and m- Let Aq and ag 
denote the positive constant as in i2S\} and Then 

(i) there exists a Lipschitz continuous function w : R™ H so that a;(a(a;)) = x for all x G H. For any 
1 < P, 9 < oo, uj has an upper Lipschitz constant Lip(w)p_g between (R"*, H-Hp) and {H, Dq) bounded by: 


Lip(w)p,g < 


•2'! 2 • max(l, m= p ) for q < 2 


’ • niax(l, m 2 p ) for q > 2 
Explicitly this means: for q <2 and for all c,d G R"*; 

Dq{ui{c), U!{d)) < ■ 2 5 ~ 2 . niax(l, m^~p) ||c — d\\^ 


+ n i_l 


whereas for q > 2 and for all c,d G 


Dq{ui{c),ui{d)) < 




3V2 + 2^+p i_i 




2 9-inax(l,m2 p)||c—d|l 


4 + 3^2 . 


In particular, for p = 2 and q = 2 its Lipschitz constant Lip(w) 2,2 bounded by 

D 2 {uj{c),uj{d)) < ||c- d|| 


(27) 


(28) 


(29) 


(30) 


(ii) there exists a Lipschitz continuous function tjj : R™ H so that 'tf{P(x)) = x for all x G H. For any 
^ F p,q < oo, if has an upper Lipschitz constant Lip('0)p^g between (R™, IHIp) and {H, dq) bounded by: 


Lip('0)p.9 < 


^±^1.2<i 2 .jnax(l,m 2 p) forq<2 
max(l, m2 “p ) for q > 2 

Explicitly this means: for q <2 and for all c,d G R"*; 

dq{ip{c),if{d)) < ^ ■ 2p“^ • max(l,m^“p) ||c- d||^ 


\/ao 


whereas for q > 2 and for all c,d G 


dq{ip(c),ip{d)) < 


3 + 2 


l + i 


,_ max(l,m 2 p)||c-d|l 


4 + 3+2 . 


In particular, for p = 2 and q = 1 its Lipschitz constant Lip('!/)) 2 ,i bounded by 

di{if{c), if{d)) < ||c - d|| 

y/af) 


(31) 


(32) 


(33) 


(34) 
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The proof of Theorem 13.31 presented in Section [31 requires construction of a special Lipschitz map. We 
believe this particular result is interesting in itself and may be used in other constructions. This construction 
is given in |12j for the case p = 2. Here we consider a general p and give a better bound for the Lipschitz 
constant. We state it as a lemma. 

Lemma 3.4 Consider the spectral decomposition of any self-adjoint operator A in Sym(iL), say A = 
^rn{k)Pk, where Xi > X 2 > ■ ■ ■ > Xn are the n eigenvalues including multiplicities, and Pi,...,Pd are 
the orthogonal projections associated to the d distinct eigenvalues. Additionally, m(l) = 1 and m{k + 1) = 
m{k) + r(k), where r(k) = rank{Pk) is the multiplicity of eigenvalue Xm{k)- Then the map 

7T:Sym{H)^ , 7 t{A) = (Xi - X2)Pi (35) 

satisfies the following two properties: 

(i) for 1 < p < 00 , it is Lipschitz continuous from (Sym(iJ), H-Hp) to H-Hp) with Lipschitz constant 

less than or equal to 3 + 2^^~p; 

(ii) 'k{A) = A for all A G 

Remark 3.5 Numerical experiments suggest the Lipschitz constant of tt is smaller than 5 for p = 00 . On 
the other hand it cannot he smaller than 2 as the following example shows. 

Example 3.6 If A = ^ ^ H ^, then 7r(H) = and tt{B) = ■ Here we have 

||7r(H) -7r(H)||^ = 2 and \\A- B\\^ = 1. Thus for this example ||7r(H) - 7r(H)||^ = 2||H- H|l^. 

It is unlikely to obtain an isometric extension in Theorem 13.31 Kirszbraun theorem |36j gives a sufficient 
condition for isometric extensions of Lipschitz maps. The theorem states that isometric extensions are 
possible when the pair of metric spaces satisfy the Kirszbraun property, or the K property: 

Definition 3.7 The Kirszbraun Property (K): Let X and Y he two metric spaces with metric d^ and dy 
respectively. {X,Y) is said to have Property (K) if for any pair of families of closed balls {B(xi,ri) : i G /}, 
{B(yi,ri) : i G I}, such that dy{yi,yj) < dx{xi,Xj) for each i,j G I, it holds that B(xi,ri) ^ 0 
(^B{yi,r,) ^ 0. 

If {X, Y) has Property (K), then by Kirszbraun’s Theorem we can extend a Lipschitz mapping defined on a 
subspace of X to a Lipschitz mapping defined on X while maintaining the Lipschitz constant. Unfortunately, 
if we consider {X, dx) = (M™, H-H) and Y = H, Property (K) does not hold for either Dp or dp. 

Property (K) does not hold for H with norm Dp. Specifically, (R™,®”/ ~) does not have Property K. 

Example 3.8 We give a counterexample for m = n = 2,p = 2: Let yi = (3, 1), ^2 = (~li 1); Va = (0,1) 
be the representatives of three points yi, y 2 , ys in R^/ ~. Then D 2 {yi,y 2 ) = 2-\/2, ^ 2 ( 2/25 2 / 3 ) = 1 and 
T> 2 {yi,y 3 ) = 3. Consider xi = (0,0), X 2 = (0,— 2-\/2), x^ = (—l,—2y/2) in R^ with the Euclidean distance, 
then we have \\xi — a; 2 || = 2-\/2, ||x 2 — =1 and ||xi — xaH = 3. For ri = Vd, r 2 = 2 — \f2, ra = \f%— \f2>, 

we see that (1 - y/2, 1 + \/2) G flLi B{xi,rf) but B{yi,ri) = 0. To see HLi B{yi,ri) = 0, it suffices to 
look at the upper half plane in R^. If we look at the upper half plane H, then B(yi,ri) becomes the union 
of two parts, namely B{yi,ri) U H and B{—yi,ri) U H, and B{yi,ri) becomes B{yi,ri) for i = 2, 3. But 
{B{yi,ri) U H)r\ 5(1/2, ^’ 2 ) = 0 and (5(-yi,ri) U H) H 5(y3,r3) = 0. So we obtain that B{yi,ri) = 0. 
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Property (K) does not hold for H with norm dp. Specifically, (R™,C"/ ~) does not have Property K. 
The following example is given in m- 

Example 3.9 (ttH) Let m be any positive integer and n = 2, p = 2. We want to show that {X,Y) = 
(]R'",C"'/ does not have Property (K). Let yi = (1,0) and y 2 = (0,-\/3) be representitives of yi, y 2 € Y, 
respectively. Then d\{yi,y 2 ) = 4. Pick any two points Xi, X 2 in X with ||a;i — X 2 \\ = 4. Then B{xi,2) and 

B(x 2,2,) intersect at x^ = {xi +X2)/2. £ X. It suffices to show that the closed balls B{yi,2) and B{y2,2) 

have no intersection in H. Assume on the contrary that the two balls intersect at j/ 3 , then pick a representive 
of Vs, say i /3 = (a, b) where a, b € C. It can be computed that 

di{yi,yz) = |a|^ + - 2 |a|^ + 2 \b\^ + 2 |a|^ + 1 (36) 

and 

<^ 1 ( 2/27 2 / 3 ) = |a| +1^1 +6|a| — 6 |6| + 2 |a| |6| +9 (37) 

Set di[yi,yfj = di{y 2 ,y 3 ) = 2. Take the difference of the right hand side of L36\) and fSTj ), we have 
l&l —|a| =landthus\b\ >1. However, the right hand side of ^3^ can be rewritten as {\a\ +|6| —1)^+4 |5| , 
so di{yi,y 3 ) = 2 would imply that \b\^ < 1/2. This is a contradiction. 

Remark 3.10 Using nonlinear functional analysis language f\14V Lemma \S.4 \ can be restated by saying that 
is a 5-Lipschitz retract in Sym(iJ). 

Remark 3.11 The Lipschitz inversion results of Theorem \3.3\ can easily be extended to systems of quadratic 
equations, not necessarily of rank-1 matrices from the phase retrieval model considered in this paper. 

4 Proofs of results 

4.1 Proof of results in Section [2] 

We start by proving Theorem 12.21 
Proof of Theorem 12.21 

(i) First we prove the following lemma. 

Lemma 4.1 Fix x £ C" and z £ C". Let f = j(a;) and f = j(z) be their realifications, respectively. 
Let ^0 G ^ {j(i) S : X € x} be a point in the equivalency class that satisfies D 2 (x, z) = ||^o ~ Cll- 

Then it is necessary that 

(eo,JC) = 0 (38) 

and 

(fo,C)>0 (39) 

where J is defined as in m- 
Proof: For 0 £ [0, 27r) define 

U{9) := cos(0)/ + sin(0) J . (40) 

Then it is easy to compute that 

iiU^x) = um • (41) 
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Therefore, 


D 2 {x,z)= min \\U(9)^-C\\ = ||C|| + ||CII “ 2 max (;7(6»)e, C) ■ 

ee [0,271-) ee[o,2'n-) 

If {U{9)^, () is constantly zero, then we are done. Otherwise, note that 

inax (C/(d)^,C) = + 

u ^[0,2 t 7) 


and the maximum is achieved at 0 = do if and only if 


and 


Now we can compute 


cos(do) 


(e,C) 


sin (do) 


{JtO 

{XX)^ + {XX)¥ 


XoXO 


{U{9oX,JC) 

cos(do)(^, JC) + sin(do)(JC, JC) 


(^,0 


(^,0 




(J^X) 


ixxr + ij^xyr 


-{-XX) + 


ixxy + {Jixyy- 

(Jix) 


-{J^XO 


ixxx + ixxyy 


-XX) 


= 0 . 


So we get (1551) . (1551) is obvious. 

Q.E.D. 

Now we come back to the proof of the theorem. Denote 


p{x,y) 


||a(3;) - a{y)f 
D2 (x,2/)2 


x,y gC^, x^ y. 


We can represent this quotient in terms of ^ and 77 . It is easy to compute that 


p{x,y) = PX,y) 


T,T=i{^XX) + {^kV,v) - 2\/($fc^,0($fc7?,7?) 

Uf + \\vf-2VX,v)^ + XXv)^ 


(42) 


(43) 


(44) 


(45) 


(46) 


Fix r > 0. Take that satisfy D 2 {x,z) = ||^ — Cll < ^ and D 2 {y,z) = |1?7 — Cll < Let 

/i = (^ + 77)72 and V = X — 77 ) 72 . Then ||z/|| < r. Note that for r small enough we have that ||/i|| > ||z/|| 
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and that 0 ^ 0- Thus 


^ Sr=i(^fc(M + + t^) + ($fc(M -v).[^-v)- 2\/($fc(M + pt + t^)(^fc(A^ -v),[i-J) 

11/^ + + \\^^ - ^\\^ - 2\/(Ai + M - I'Y + {^i, + v,J{^i- v)y 

^ Sr=l(^fcM> ^^) + ^) - \/M) + - 4($fc/i, Z/)^ 

ll^f + IIHI^ - \/||/.r + |kf-2||Mf llz/f+ 4 (m,^^)^ 

^ Sfc:$,C#o(4’fc/i, A^) + v) - ^/^{^k^J-, m) + (4>fct/, t/))^ - 4(^fc^, vY 

llA^f + IIHI^-VM" + IIHl"-2||Mf Ikf 

Efc:$,c#o(4'fcM, A^) + m) + (4»fct^, 




ii 



Wki^^ 




(4>fcAi,Ai) ”*" 


. {‘S>klJ.,i^y 
WkV^ 



= ^(5(/iMz.)+0(||z.f) . 


Note that 


|(J/i,j^)| = |(J/i,J^) - (TC,i^)l < I|Tm- TCII ||J^|| = llti- Cll 11*^11 


(47) 


since {JQ,v) = 0 by Lemma ITTI Also, ||^ — C|| < r. Therefore, 


\\PjA\ 


I(Tm, v)\ ^ \{J^J-.v)\ ^ r\\v\\ 


WJt^W M - IIa^II 


and thus 
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As a consequence, we have 


P{^,V) = 


> 


{S{lJL)Pj- V,Pj~ v) y 

0(lk|| ) 

||z.f 

{S{^,)P^^u,P^^v) ' 


1 - 




+ 0{r^) 


> 1 - 


^2n-l (‘5(/x)) + 0(r^) . 


Take r —> 0, by the continuity of eigenvalues with respect to matrix entries we have that 

A(z)>A2„-i(5(C)) . (48) 

On the other hand, take £’ 2 ™-! to be the unit-norm eigenvector correspondent to A2n-i(‘5(C))- For 
each r > 0, take ^ = ( + §£ 2^-1 and rj = ( — §£ 2 n-i- Then 

p{x,y) = £(C,?7) = A2n-i(5(C)) • 


Hence 


Together with (H5)) we have 


A(z) < X 2 n-l{S{C)) ■ 


A(z)=A2„-i(5(C)) . 

(ii) Assume on the contrary that Aq = 0, then for any iV G N, there exist xn, Un ^ H ior which 

, ||a(a;iv) - cx{yN)\f ^ 1 

P(«,w)= ^N- 


(49) 

(50) 

(51) 


Without loss of generality we assume that ||a:jv|| > ||yAr|| for each N, for otherwise we can just swap the 
role of a:AT and y^- Also due to homogeneity we assume ||a:Ar|| = 1. By compactness of the closed ball 
£ 1 ( 0 ) = {x £ H : ||a;|| < 1} in £ = C", there exist convergent subsequences of {xn}^^^ and {yN}N£N, 
which to avoid overuse of notations we still denote as {xn}ngn Xq £ H and {yN}NGN yo £ H. 

Since ||a;o|| = 1 we have from (i) that A{xo) > 0. Note that D 2 {xN,yN) < + ||yAr|| < 2, so by 

m we have ||a(a;Ar) — a(?/Ar)|| —>■ 0. That is, ||Q;(a;o) ~ Q^(yo)ll = 0. By injectivity we have Xg = yo in 
H. By ProDOsition l2.2f i'). 

p{xN,yN) > A(a;o) - > l/N 

for N large enough. This is a contradiction with (ICTl . 

(iii) The case z = 0 is an easy computation. We now present the proof for z ^ Q. First we consider 
p{x, z) = P{^, () as defined in (|4^ . Fix r > 0. Take ^ G that satisfy D 2 {x, z) = ||^ — C|| < r. Let 

d = X — z and <5 = j(d) = ^ — C- Note that 


lief+ Iicf- 2 v'(e,C)^ + (e,JC)^ 
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The numerator is equal to 


C) + 2($fcC, <5) + 5) + ($feC, C) - 2v^(($fcC, C) + 2($fcC, <5)) • i^kC, 0 

fc=l 


= ^ 2($fcC,C)+2($fcC,<5) + ($fc(5,<5)+2(<i>fcC,C)[l + 

fe:$fcC#0 


{<^kCXV 


1 4($fcCC)^(^fcC<^)^ 
8’ ($feC,C)^ 


k:^kC =0 


= E 

k-.-s-kC^o 


{<^kC 6 )^ 

i^kCO 


+ E {^k6,S) + 0{\\5f); 

k:^kC =0 


the denominator is equal to 


2||Cf+ ||5f+2(C,5) 
= ll<5f+ 0(||<5f) 


, IlCf (C,5) + ^(C,<5) + i(JC,<5)^ 
2IICII (^1+ 


4||Cf 

8||Cf 



in which we used Lemma [4.II to get {JC,^) = 0. 


Take r —>• 0, we see that 


A{z) > \ 2 n -1 + 


k-.{z,fk)=0 




(52) 


Let E2n-i be the unit-norm eigenvector correspondent to \2n-1 {S{C) +J2k-{z,fk)=o^^)- Note that 

{JC: E2n-i) = 0 since 5(C) JC = 0 and = J^k( = 0 for each k with {z,fk) = 0. Take C = 
C + i^ 2 n-i for each r, we again also have 


Therefore 


A{z) < X2n-1 I 5(C) + d’fe 

k-{zjk)=0 


A{Z) = X2n-1 I 5(C) + E 

fc:(z./fc)=0 


(53) 


(54) 


(iv) Take z = 0 in (iii). 

(v) B(z) can be computed in a similar way as in (iii) (in particular, the expansion for T’(C, C) is exactly 
the same). We compute B{z). 5(0) is computed in [T3], Lemma 16. Now we consider z ^0. Use the 
same notations as in (ITO . Fix r > 0. Again, take rj G that satisfy D 2 {x, z) = ||C — Cll < »' and 
D 2 {y,z) = \\r] — Cll < r. Let fJ- = {^ + ? 7)/2 and v = {i — ri)/‘^- Also let (5i = C — C and 82 = V — C- 
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Recall that 


P{^,V) = 


Er=i(^feC>o+ {^kv,v) - ‘2^/Wk^x)Wk^^ 


Uf + Wvf-^VW^WT^JW 


= E 


(5>fc^, 0 + i^kV, V) - 2\/($fc^,0(®fc^>?7) 


ti II^II + ihir - 2.MWTWJW 

Now we compute it as J2T=i = 'Ek-.<s>^c^o + Efe:$,c=0' 

0 + {^kV, v) - 2\/($fc^, 0{<^kV, v) 

k:i:t^o ii^f + M" - 2v/(e,^)" + (e,j^)^ 


= E 

fe:$fcC#0 


^l) + I/) - ^y (($fcM, //) + :/))2 - Z/)2 


IVIIHI^- 


(55) 


The computation for its numerator is the same as in (i). We get that the numerator is equal to 

2(5(m)z.,./)+0(||H|") . 

Since /r 7 ^ 0, the denominator is equal to 






1 - 


iiHr , 2(/x,ji.)2 


+ 


+o(iiHr) 


=2 M- 


2 2 {Jfi,i^y 


+ o(iiHr) 


= 2 \Wf + 0 {\Wf) bygll). 

Also we can compute using the denominator as above [note that i' = {Si — S 2 )/ 2 ] that 


E 

k-.^kC =0 


i^kto + {^kV,v) - 2v^($fcg,^)($fc7?,?7) 

mf + \\vf-2iMw+{(JW 


= E 

fc:$fcC=0 






U ll«i-&ll" + o(M‘) 


Now put together (1551) . (1551) and (1571) . we get 


P{^,V) = 


{S{n),.,^) + 0{\Wf) , ^ (lhE^i||-|hE<52||) 


¥\? + o{\\vt) 


E 

fc:'I’fcC=0 


||5i-52f+ o(||Hl") 


Note that 






^ — (^ 2 ),(5i — < 52 ) 


(56) 


(57) 


(58) 
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since it is equivalent to 


(59) 


($fel 5 i, 5 l)($fclJ 2 , <^2) > ^2))^ 

which is the Cauchy-Schwarz inequality. Therefore we have that 

^ +0(|kf) 

< - „ „2 „4, - < I ‘5(/r) 


wr+omn 


$fc|+0(r2). (60) 

fe:‘I'fcC=0 


Take r —?► 0 we have that 


B{z) < Ai 5(0 + E 

\ fe:«'fcC=0 

Again we get the other direction of the above inequality by taking ^ = ( + ^Ei and r] = ( — ^Ei for 
each r > 0 where Ei is the unit-norm eigenvector correspondent to Ai ^5(0 + J2k-{z fk)=o Note 
that for each r, the equality in (l5^ holds for this pair of ^ and 77 . 

(vi) Take z = 0 in (v). 

Now we prove Theorem l2.6l 

Proof of Theorem 12.61 

Only the first two parts are nontrivial. We prove them as follows. 

Fix z G C”. Take x = z + di and y = z + d 2 with ||di|| < r and ||d 2 || < r for r small. Let u = x + y = 
2z + di+ d 2 and v = x — y = di — d 2 . Let fj, = 2( + 61 + S 2 G and 1 / = Si — S 2 & be the realification 
of u and v, respectively. Define 


p{x,y) = 


\\l3{x)-PiyW 


di{x,yY 

By the same computation as in [7], Section 4.1, we get 

{TZ{2( -I- (5i -I- (52)(<5 i — ^ 2 ), <5i — S 2 ) 


p{x,y) = Q(C;di,(52) := 


|| 2 C -b di -I- ( 52 II {Pj~( 2 C+Si+S 2 )^^'^ ~ '^2),'^! - S2) 
Since J{2C, -b di -b ^2) G ker 'R-{2(^ -b di -b ^2), we have 

Q(C-§ § ) — + ^1 + '^2)Pj‘(2,^+j^-|_fe)(<^l ~ ^ 2 ), Pj(^2C+Si+S2)^^^ ~ 

’ P<: + Sl+S 2 \f{Pj-^,^^s,^S 2 )(^l-S 2 ), 6 i- 62 ) 

Now let (5 = (5i -b (52 and i/ = 61 — 62 - Note the set inclusion relation 


(61) 


(62) 


(63) 


^2 s 
c { 5 i ,52 g : 
c {di,d 2 e: 


j2n 

: 

D2n , 


Thus we have 


„inf (3(C;5 i,<52) < 

||5||<2r 

\\v\\<2r 

vXJ{2C,+&) 


||d||<^, ||H|<^, zzT J(2C + <5)} 

||(5i||<r, ||(52|| < r, 1 /T J(2C-b(5)} 
||(5|| < 2r, ||zz|| <2r-, z/T J(2C-b^)} 


inf Q(C;di,(52) < „^inf Q(C;(5i,(52) 

||5i||<r \\S\\<r/2 

|F 2 ||<r M\<r/2 

i2±J(2C+S) i2±J(2C+S) 
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That is, 


^2n-l{T^{2C, + S)) 

5ir<2r ||2C + (5f 


< inf QiC; 61 , 62 ) < inf ^ 

ll^illO' ll'5||<r/2 |12C + ^|| 

||(52||<r 
i/_LJ(2C+<5) 


Take r —>■ 0, by the continuity of eigenvalues with respect to the matrix entries, we have 


That is. 

Now consider 


A2„_l(7^(C))/ IlCf < a{z) < A2„_l(7^(C))/ | 


a(z)=A2„-l(7^(C))/||Clr ■ 


pix,z) = 


For simplicity write <5 = (5i. We can compute that 

(7^(2C + <5)<5,5) 


di{x,z)'^ 

_ (7e(2C + 5)Pj(2c+,) ^’Pjl2(:+s)^) 


Note that 


p(x, z) = Q(C;6) = „ , — „ , 

l|2C + {Pj-^2c+s)^, 6) ||2C + <5f <5) 


inf Q{C, 6 ) > inf inf Q{C, 6 ) = inf A2n-i(P(2C + 5)) 

||<5||<r IklKr ||i5||<r lk||<r 

<5_LJ(2C+5) <5_LJ(2C+5) 


(64) 

(65) 

( 66 ) 


Take r —>■ 0 we have that 


d(z) > A2 „-i(P( 2C))/ ||2Cf = A 2 „-i(P(C))/ IlCf 


On the other hand, take e 2 n-i to be a unit-norm eigenvector correspondent to A2n-i(7?.(2C)). Then by 
the continuity of eigenvalues with respect to the matrix entries, for any e > 0, there exists t > 0 so that 
6 = te 2 n-i satisfy 

(7Z(2( + 6)6,6) ^ (^(20) + e (67) 




and from there we have 


Therefore, 


d(z) < A 2 „-i(P( 2 C))/ ||2Cf = A2u-i(7Z(C))/ IlCf 


a(z) = A2n-i(7Z(C))/IICir 


( 68 ) 


(69) 


In a similar way (replacing infimum by supremum) we also get b(z) and b(z) as stated in the theorem. 

Q.E.D. 
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4.2 Proof of results in Section [3] 

We start by proving Proposition 13.II 

Proof of Proposition [37T] 

(i) For Dp obviously we have Dp[x,y) > 0 for any x, y G H and Dp{x,y) = 0 if and only ii x = y. We 

also have Dp{x,y) = Dp{y,x) since \\x — aj/Hp = ||?/ — for any x, y G H, |a| = 1. Moreoever, 

for any x,y, zG H, if Dp{x, y) = ||x - ayW^, Dp{y, z) = \\z - by\\, then 

Dp{x, z) < ||x - ab~'^z\\^ = \\bx- az\\p < \\bx - aby\\p + \\aby - az\\p = Dp{x, y) + Dpiij, z). 

Therefore Dp is a metric. 

dp is also a metric since H-Hp in the definition of dp is the standard Schatten p-norm of a matrix. 

(ii) For p < g, by Holder’s inequality we have for any x = (xi,X 2 , ...,x„) G H = that 1^*1^ — 

^ Thus ||x||p < lla^H^. Also since IHIp is homogeneous, if we assume ||x||p = 1 

we have < J2i=i \d\^ = 1- Thus ||x||^ < ||x||p. Therefore, we have Dq{x,y) = ||x - aiyW^ > 

II 2 ; _ aiyWj^ > n^p~p'>Dp{x, y) and Dp{x, y) = ||x - a 2 y\\p > ||x - a 2 y|l^ > T>,(x, y) for some Oi, 
02 with magnitude 1. Hence 


Dq{x,y) < Dp[x,y) <n'p p^Dq{x,y) 

We see that (Tlp)i<p<oo are equivalent. The second part follows then immediately. 

(iii) The proof is similar to (ii). Note that there are at most 2 ai’s that are nonzero, so we have 2^p~p'^ 
instead of rSp~p\ 

(iv) To prove that Dp and dq are equivalent, we need only to show that each open ball with respect to Dp 
contains an open ball with respect to dp, and vise versa. By (ii) and (iii), it is sufficient to consider 
the case when p = q = 2. 

First, we fix X G H = C”, r > 0. Let R = min(l, ( 2 ||x|| ''+i)n ' -^ Then for any y such that D 2 {x, y) < R, 
we take y such that ||x — y\\ < R, then VI < i,j < n, \xixj — ytyj] = |xi(^ — y]) + (xj — yi)yj\ < 
\x^\R + R{\xi\ +R) = i?(2|xi| +R) < i?(2|xj| + 1) < ;^. Hence d 2 {x,y) = ||xx* - yy *\\2 < = r. 

On the other hand, we fix x G 7L = C", R> 0. Let r = R? 1^/2. Then for any y such that d 2 {x, y) < r, 
we have 

{d 2 {x,y)f = ||xf + \\yf - 2|(x,y)p < ^ (70) 


But we also have 


{D 2 {x, y)f = niin ||x - apf = 

|a| = l 


{x,y) 

X - 7^ - -y 


= ||xf+ ||yf-2 |(x,2/)| 


So 


l(a;,y)r 

(Z? 2 (i,y))" = Ikf + llpf + 2 ||xf \\yf - 4(||xf + \\yf)\{x,y)\ + 4\{x,y)\^ 


(71) 


(72) 


19 








Since \{x,y)\ < ||x|| \\y\\ < (HxH^ + ||y||^)/2, we can easily check that {D 2 {x,y))'^ < 2{d2{x,y)Y < R'^- 
Hence D 2 {x,y) < R. 

Thus Z ?2 and c ?2 are indeed equivalent metrics. Therefore Dp and dq are equivalent. Also, the imbedding 
i is not Lipschitz: if we take x = (a;i, 0,..., 0) G C", then D 2 {x, 0) = |a;i|, d 2 {x, 0) = \xi\^. 

(v) First, for p = 2, ior x ^ y va H — {0}, we compute the quotient 

||Ka(a;) - Ka(y)f 

p{x,y) = -ppT—72- 


\\yf + Ikf \\yy*f - 2 Ik|| Ibll trace(a;x*TO*) 

\\x\\^\\yf + \\xf\\y\\^- 2 \\xf\\yf\x*y\ 

^ 2 Ikll ||x/|| {\\x\\ ||y|| \x*y\ - tTace{xx*yy*)) 

\\x\f\\yf + \\xf\\yf- 2 \\xf\\yf\x*y\ 

^ 2 (Ikll ||?/|| \x*y\- tTace{xx*yy*)) 

INni2/ll + INIIIyf-2|N|||2/|||x*y| 


\\Ka{x) - Ka{y)\\ 

D 2 [x,yY 

\\x\\~^xx* - \\y\\~^yy* 
\\xf + \\yf- 2 \{x,y)\ 


where we used ||xa;*|| = ||a;||^. For simplicity write a = ||x||, b = ||?/|| and t 
have a > 0, 6 > 0 and 0 < t < 1. 


Now we have 


p{x,y) = 1 + 


2 {abt — abt^) 
(j 2 _|_ ^2 _ 2 abt 


Obviously, we have p{x,y) > 1. Now we prove that p{x,y) < 2. Note that 


\{x,y)\ ■ (Ikll II2/II) ^ We 


1 + 


2 {abt — abt^) 
q _2 _|_ j ,2 _ 2 abt 


< 2 + 6^ — ^abt + 2 abt^ > 0 


But 


+ b^ — Aabt + 2 abt^ > 2 ab — Aabt + 2 abt^ = 2 ab{t — 1)^ > 0, 


so we are done. Note that take any x, y with (x, y) = 0 we would have p{x, y) = 1. On the other hand, 
taking ||x|| = ||y|| and let t —>■ 1 we see that p(x, y) = 2 — e is achievable for any small e > 0. Therefore 
the constants are optimal. The case where one of x and y is zero would not break the constraint of 
these two constants. Therefore after taking the square root we get lower Lipschitz constant 1 and 
upper Lipschitz constant ^/2. 

For other p, we use the results in (ii) and (iii) to get that the lower Lipschitz constant for is 
min(22 p,7t,p 2 j and the upper Lipschitz constant is v2 max(n 2 p,2p ^). 


(vi) This follows directly from the construction of the map. 


(vii) This follows directly from (v) and (vi). 

Q.E.D. 
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Next we prove Lemma [3.41 
Proof of Lemma 13.41 

(ii) follows directly from the expression of tt. We prove (i) below. 

Let A, B G Sym(iJ) where A = X]fe=i ^m{k)Pk and B = X]fe'=i P-m(k')Qk'- We now show that 


\7riA)-7riB)\\<{3 + 2^H)\\A-B\\ 


(73) 


Assume Ai — A 2 < fJ-i — lJ- 2 - Otherwise switch the notations for A and B. If /ri — = 0 then 7 r(A) = 

Tr{B) = 0 and the inequality (17^ is satisfied. Assume now m — fi 2 > 0. Thus Qi is of rank 1 and therefore 
llQilIp = 1 for all p. First note that 


7r(A) — tt{B) — (Ai — A2 )Pi — — ^2)Qi — (Ai — A2)(Ti — Qi) + (Ai — — (A 2 — P 2 ))Qi 


(74) 


Here HPiUgo = llQilloo = Therefore we have ||Pi — ^ 1 since Pi, Qi > 0. From that we have 

\\Pi-Qi\\p<2-^. 

Also, by Weyl’s inequality we have |Ai — pi\ < ||A — for each i. Apply this to i = 1, 2 we get 
|Ai — Pi — (A 2 — ^ 2)1 < |Ai — ^i| + IA 2 — P 2 \ < 2 ||A — B\\^. Thus |Ai — ^i| + IA 2 — P 2 \ < 2 ||A — B\\^ < 
2\\A-B\\^. 

Let g :=Ai — A 2 , S := ||A — P|| then apply the above inequality to (ITU) we get 


| 7 r(A) - 7 r(P)|| < g ||Pi - Qi\\ +2S <2pg + 2S 


(75) 


If 0 < 5 < (2 + 2 p)S, then || 7 r(A) — tt{B)\\^ < {2^'^p + 3)5 and we are done. 

Now we consider the case where g > (2 + 2~p)5. Note that in this case we have 5 < g/2. Thus we have 
|Ai — pi\ < 5/2 and IA 2 — P 2 \ < g/2. That means pi > (Ai + A 2)/2 and p 2 < (Ai + A 2 )/ 2 . Therefore, we 
can use holomorphic functional calculus and put 


Pi =- — (f 

2Tn X 


RAdz 


and 


Qi = - ^ f Rsdz 


(76) 


(77) 


where Ra = {A — zl)~^, Rb = {B — zl)~^, and 7 = ^{t) is the contour given in the picture below (note 
that 7 encloses pi but not P 2 ) and used also by [37]. 

Therefore we have 

^ " .- (78) 

11 


\Pi - Qillp I \\{Ra - RB)hit))l wmt 


Now we have 

(Ra - Rb){z) = Ra{z) - (/ + Ra{z){B - A))-^Ra{z) = {Ra{z){B - A)TRa{z) (79) 


n>l 


since for large L we have ||Pa(^:)(P - A)||^ < ||P.4(^:)IL 11^ “ ^llp < dist(z,p(A)) ^ T < 
p{A) denotes the spectrum of A. 
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Therefore we have 

\\{ra - RsKimip < E w^AiimC^ ii^ - ^^ii; 


n>l 


\\RAhmi\\A-B\\ 


< 


\\A-B\\ 


(80) 


1 - ||i?A(7(0)lloo P - -Slip p{A)) 


(2 


i+i 


1 ) 


since dist{'y{t), p{A)) > g/2 for each t for large L. Here we used the fact that if we order the singular values 
of any matrix X such that cri(^) > cr 2 (^) > • • •, then for any i we have ai{XY) < ai{X)ai(Y), and thus 
for two operators X,Y£ Sym(i?), we have |lXy||p < ll^llp- 

Hence by (1751) and (1501) we have 


17^1 - Qill < (2^+2-!)^^—^ 


-W{t)\dt 


rr dist‘^{j{t),p{A)) 

By evaluating the integral and letting L approach infinity for the contour, we have as in m 


dist‘^{j{t), p{A)) 


\l'{t)\dt = 2 


rdt = 


Hence 


0 + 
_J\A-B\\^ 2n 


'4 

^2t\ 

— arctan 

- 
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9 


\Pi - QiL < (2" + 2-')2-^ . — = (2i+p + 1) 


(81) 

(82) 

(83) 


Thus by the first inequality in (TTSl) and (1551) we have |l7r(H) — 7r(i?)||p < (3 + 2^^p)5. 

We have proved that ||7r(H) — 7r(i?)|| < (3 + ||H — B\\ That is to say, tt : (Sym(i7), H-H ) 


{S^’^{H), IHlp) is Lipschitz continuous with Lipschitz constant less than or equal to 3 - 

Q.E.D. 


H+f; 


Now we are ready to prove Theorem 13.31 
Proof of Theorem 13.31 

The proof for a and /3 are the same in essence. For simplicity we do it for /3 first. 
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We construct a map V' : (R™, IHIp) —>■ {H,dq) so that ^(/3(x)) = x for all x € H, and tp is Lipschitz 
continuous. We prove the Lipschitz bound (l27l) which implies (l26l) for p = 2 and g = 1. 

Set M = I3{H) C M™. By hypothesis, there is a map ipi : M ^ H that is Lipschitz continuous and 
satisfies ipi{f]{x)) = x for all x £ H. Additionally, the Lipschitz bound between (M, 11 - 112 ) (that is, M with 
Euclidian distance) and {H,di) is given by 

First we change the metric on H from di to d 2 and embed isometrically H into Sym(iL) with Frobenius 
norm (i.e. Euclidian metric): 


(M,|H| 2 ) ^ (Lf,di) ^ {H,d 2 ) ^ (Sym(iL),|H|^,) 


(84) 


where * 1 , 2 ( 2 :) = x is the identity of H and np is the isometry (1^ . We obtain a map 1^2 ■ {M, 11 - 112 ) — 
(Sym(iL), of Lipschitz constant 


Lip(*/’2) < Lip(i/(i)Lip(*i,2)Lip(H;^) 


1 

’ 


where we used Lip(*i, 2 ) = Lf 2 n = 1 by (ITOl) . 

Kirszbraun Theorem [55] extends isometrically 1^2 from M to the entire R'” with Euclidian metric ||-||. 
Thus we obtain a Lipschitz map i /:2 : (R'",||-||) —>■ (Sym(7L), ||-||^^) of Lipschitz constant Lip (^ 2 ) = Lip(i/: 2 ) < 
so that 'ip2{l3{x)) = [ 2 :, a:] for all x £ H. 

The third step is to piece together ip 2 with norm changing identities. 

For q < 2 we consider the following maps: 


(R™,||-||p) ^ (R’”,||-|| 2 ) A (Sym(iJ),||-||^J ^ ^ (^.'^ 2 ) ^ {H,dq) (85) 

where jp ,2 and * 2,5 are identity maps on the respective spaces that change the metric. The map ip claimed 
by Theorem 13.31 is obtained by composing: 

^ (H^dq) , Ip = i2,q ■ ■ tp2 ■ jp,2 

Its Lipschitz constant is bounded by 

Lip(i/i)p,g < Lip(jp, 2 )Lip('!/: 2 )Lip( 7 r)Lip(K 7 ^)Lip(* 2 ,q) < max(l, m^~p)—— ■ (3 + 2^2) ■ 1 - 29 “ = 


Hence we obtained ((32ll . The other equation (I26p follows for p = 2 and g = 1. 

For g > 2 we use: 

(R™,||-||^) ^ (R'",|H| 2 ) ^ (Sym(iL),|M|^,) ^ (Sym(iJ), ||-||,) ^ (^^’“(iL), |M|,) ^ iH,dq) (86) 

where jp ,2 and / 2 ,g are identity maps on the respective spaces that change the metric. The map ip claimed 
by Theorem 13.31 is obtained by composing: 

i/::(R’",||-||p)^(iL,dg) , iP = n-^^-7r-l2,q-iP2-Jp,2 

Its Lipschitz constant is bounded by 

Lip{ip)p^q < Lip(jp,2)Lip('02)Lip(/2,g)Lip(-7r)Lip(«:o^) < max(l,)—^ . 1 • (3+ 2 ^+ 9 ) . 1 

y/a-o 
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Hence we obtained (1551) . 

Replace by a, '0 by and np by in the proof above, using the Lipschitz constants for in 
Proposition 13.11 we obtain (1551) and (1551) . 

Q.E.D. 
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